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Abstract. In three previous papers (Pelat 1997, 1998 and Moultaka & Pelat 2000), we set out an inverse stellar population 
synthesis method which uses a database of stellar spectra. Unlike other methods, this one provides a full knowledge of all 
possible solutions as well as a good estimation of their stability; moreover, it provides the unique approximate solution, when 
the problem is overdetermined, using a rigorous minimization procedure. In Boisson et al. (2000), this method has been applied 
to 10 active and 2 normal galaxies. 

In this paper we analyse the results of the method after constraining the solutions. Adding a priori physical conditions on the 
, solutions constitutes a good way to regularize the synthesis problem. As an illustration we introduce physical constraints on 

the relative number of stars taking into account our present knowledge of the initial mass function in galaxies. In order to avoid 
biases on the solutions due to such constraints, we use constraints involving only inequalities between the number of stars, after 
dividing the H-R diagram into various groups of stellar masses. 

We discuss the results for a well-known globular cluster of the galaxy M31 and discuss some of the galaxies studied in Boisson 
et al. (2000). We find that, given the spectral resolution and the spectral domain, the method is very stable according to such 
constraints (i.e. the constrained solutions are almost the same as the unconstrained one). However, an additional information 
can be derived about the evolutionary stage of the last burst of star formation, but the precise age of this particular burst seems 
to be questionable. 
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1. Introduction 



■ The search of the stellar populations inside unresolved galaxies has been the aim of several studies since the sixtees. Two 
' different approaches have been adopted for this purpose: the direct approach of which methods are usually called "the evolutive 
synthesis methods" (e.g. Tinsley 1972, Chariot & Bruzual 1991, Bruzual & Chariot 1993, Leitherer et al. 1999, Fioc & 
Rocca-Volmerange 1997, Vazdekis, 1999, Bruzual & Chariot 2003) and the inverse one of which these are called "the synthesis 
methods" (e.g. Faber 1972, O'Connell 1976, Joly 1974, Bica 1988, Schmidt et al. 1989, Silva 1991, Pelat 1997,1998, Moultaka 
& Pelat 2000). 

In the first approach, one decides an a priori model for the history of the star formation occuring inside the studied galaxy, and 
by means of theoretical stellar evolutionary tracks and of a stellar database, derives quantities that are directly compared to the 
observed ones. From different input models, one retains the model that best fits the observed quantities. 



In the inverse approach, usually, no a priori model is necessary to derive the stellar populations and one uses exclusively 
the observables in order to deduce the stellar spectral types and luminosity classes by means of a minimization procedure. 
This minimization task is not a very simple one since the absolute minimum is usually difficult to find because the "objective 
function" (which is the function that one has to minimize) can rarely be minimized analytically. 

Whatever the approach is, the problem of stellar population synthesis often suffers from the lack of true solutions or, on 
the contrary, from their degeneracy (i.e. multiple solutions) and/or finally from their instability (i.e. small errors around the 
observations can induce discontinuities in the solutions). These three inconveniences have been controlled in the inverse method 
described in Pelat (1997, 1998, hereafter Paper I and II) where all the solutions are well identified and the minimization procedure 
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is rigorously treated, as well as in Moultaka & Pelat (2000) where the stability of the various solutions is analysed (hereafter 
Paper III). 

In this paper, we study the influence upon the different solutions of astrophysical constraints included a priori when searching 
for a solution. The concept of constraining the solutions in the inverse methods, has been adopted by O'connell (1976), Pickles 
(1985) and Silva (1991) but no estimation of the induced bias has been given by these authors to our knowledge. 

In the next section, we recall briefly the inverse method and its error analysis; in the third section, we describe constrained 
models. Finally, in section^ we show and discuss constrained versus unconstrained results for the globular cluster G170 located 
in M31, the LINER NGC4278, the starburst NGC3310 and the Seyfert 2 galaxy NGC21 10. In the last section we make a general 
description of the behaviour of constrained solutions in the 27 central regions of the twelve galaxies studied in Boisson et al. 
(2000), hereafter Paper IV. 



2. Description of the inverse method 

As described in Pelat (1997,1998), the present inverse method uses the equivalent widths of galactic spectra absorption lines 
as observables to be fitted by a combination of the continua fluxes and equivalent widths of a stellar lines database. The basic 
equation providing the synthetic equivalent widths is the following: 

Wsynj = „, for j = l,...,n A (1) 

where Wsynj is the synthetic equivalent width of a line at wavelength Aj, W), and are respectively, the equivalent widths 
and continua fluxes (normalized at a reference wavelength Aq) of the same lines measured in stars of class z; is the contribution 
of star class ; to luminosity at the reference wavelength Aq; n+ and nx are respectively the total number of stars considered in the 
database and the total number of lines measured in the spectra. 

To this equation, we add two physical conditions, that all the stellar contributions to luminosity kj are positive and their sum 
equals one: 

k > (2) 
J> = 1 (3) 

!=1 

Thus, having the observed set of galactic equivalent widths W ^, Sj , one searchs for the stellar contributions k, satisfying the 
previous physical conditions and minimizing the following objective function which is the square of what is called the synthetic 
distance D. The synthetic distance represents a kind of y 1 considering the W b S j as data. In fact D 2 would be exactly a^ 2 if Ft 1 
were choosen as the W b s / variances (i.e. Pj = cr^ , see Sect. 5 of Paper I for a discussion on cr WalisJ value) : 

"A 

D 2 = £(W obsj " WsynjfPj, Pj > (4) 

In this definition, Pj is a weight characterizing the quality of the equivalent widths measurements. One can eliminate this weight 
by applying the change of variables Wj <-> Wj yPj. In the following, we will consider this change of variables already made. As 
stated, the problem can either be overdetermined (i.e. there are more equations than unknowns) or underdetermined (i.e there are 
more unknowns than equations). In the first case, one finds at most a unique solution and, because of observational errors, there 
is usually not any exact solution; then the adopted solution is the approximate one which minimizes the synthetic distance. In 
Paper I, it is shown that this minimum is unique and near the true one; in addition, it is demonstrated that when signal to noise 
ratio goes to infinity, the unique minimum is exactly the true one (see Paper I). 

In the underdetermined case, one gets an infinite number of solutions which are convex combinations of particular solutions 
called the extreme solutions (see Paper II). 



Finally, the search of the error regions around the solutions has been made in Paper III, giving thus a relevance for the 
solutions and a criterium to sort the various solutions in the underdetermined case by order of merit. An "ideal" database (i.e. 
a database with an infinite spectral resolution, the largest wavelength domain and adapted for the velocity dispersion observed 
in the studied galaxy) is not degenerate. Then the analysis made in Paper III gives equivalently a condition to the database 
(depending on the quality of the observations) allowing to get a well-defined solution. Indeed, depending on the quality of the 
observations (i.e. on the size of the galactic error region), the stellar database may become degenerate. As the S/N ratios of the 
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stellar spectra are usually higher than those of the galaxy, more than one star may lay inside the galactic error region in the 
equivalent widths vector space. This situation leads to very badly determined solutions, because the different stars can not all be 
distinguished in this case. One has to eliminate such correlated stars from the database in order to obtain a well-defined solution. 
In the Appendix, the computation of the standard deviation around the synthetic distance is described. 



3. The constraints model: 

3.1. The stellar database: 

The database used comes from three different stellar libraries Serote Roos et al. (1996), Silva & Cornell (1992) and Fluks et al. 
(1994). The spectral resolution of the resulting stellar database is of 11A and the spectral domain goes from 5000A to 8800A. 
Given the galaxy velocity dispersion, the velocity broadening of the lines is equivalent to the spectral resolution of the stellar 
database. Thus no correction for velocity dispersion has to be applied. In this domain and at this spectral resolution, we have 
selected 47 line features and measured their equivalent widths (see Paper IV), some wavelength intervals have not been included 
in the database because of atmospheric absorption bands not corrected in Silva & Cornell library. The table of the line intervals 
is shown in Boisson et al. (2000). Stars included in the database have been chosen so that the H-R diagram is best represented in 
spectral types and luminosity classes for two metallicities (solar and about twice solar). 

3.2. The model: 

In general, the synthesis problem is, what is called in mathematical terms, an "Ill-posed" problem. This means that the problem 
may have no solution or a large number of solutions and/or that the solution is not stable against small deviations around the 
observation. The usual way to overcome this difficulty is to regularize the problem. This can be done by searching the solution of 
maximal entropy for example or using any other reasonable criterium. We suggest here to regularize the problem using physical 
criteria in order to obtain a unique and stable solution. This regularization procedure has already been adopted in the previous 
papers when the positivity conditions (k, > 0) were considered. In the present paper, we introduce more physical conditions in 
order to constrain the solution. 

Technically, constraining the solutions of an inverse problem comes to reduce the volume of the simplex of solutions or the 
synthetic surface as defined in Paper I (which is the set of exact solutions). The problem of such a procedure, is that the solution 
may be biased by the constraints model while, according to the philosophy of the inverse approach, it should not be altered by 
any a priori overconstrained model (otherwise it will reflect the model itself and will not provide additional information on the 
real stellar population). This inconvenience can be checked if the solutions are found to be on the border of the simplex (where 
this one has been reduced), because in this case, this would mean that the solution has been strongly constrained in order to lay 
in the reduced simplex (otherwise, it would have layed inside or on the other borders of the simplex). 
Then we could summarize the process for the search of a solution with the following scheme: 

Having the data (the spectra), one uses a model in order to derive the stellar contributions to luminosity this is the inverse 
approach. This process could, as shown in Paper II, provide an infinite number of exact solutions (in the underdetermined case) 
or an approximate solution obtained by a minimization procedure (in the overdetermined case). The obtained solutions are 
"mathematical" solutions which solve rigourously the mathematical problem. They could be stable or unstable against small 
deviations around the observation. In order to reduce the number of solutions and to insure that the solutions satisfy the physical 
conditions of the studied object, it is necessary to constrain the problem. 

Thus, the synthesis problem may be stated as minimizing the synthetic distance of equation @ also written in the following 
form: 



(where the change of variable Wj <-> Wj yPj has been applied). The set of stellar contributions k = {k\, k2, ...£«„) is submitted 
to the conditions: 



In the previous inequalities, / and u are respectively, lower and upper constant limits and the first inequality takes into account 
that all contributions are positive and less than one because of condition 0, therefore, the first n* components of vectors / and 




(5) 




(6) 
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u are, respectively, equal to zero and one. C is the constraints matrix of columns in which the first line is a vector of com- 
ponents equal to one in order to express constraint number then the component number + 1 of vectors / and u is equal to one. 

Since the expression of the synthetic distance of equation (JSJl is not linear in k, we decided to minimize instead the linear 
function ||A£|| 2 submitted to the above constraints (where A is a matrix defined in Paper I of whose components are A,-; = 
(Wobsj ~ Wji)Iji). This function happens to be the "principal synthetic distance" defined in Paper I as: 

D prin = I>obs, " Wsyn/^yn,, / S yn; = £ Ij& (7) 

j=i '=i 

As it is stated in Paper I, the problem possesses then a unique solution. Thus, using the constrained least square method, we 
can derive the unique set of stellar contributions minimizing the principal synthetic distance and satisfying the model constraints 
JSJ. Once the principal solution is at hand, we can find the "principal geometrical" one which is the nearest solution to the 
constrained principal solution and minimizes the initial synthetic distance of equation l|5}- For this purpose, we use the same 
iterative method as in Paper I where the iteration is made on the synthetic surface. At each step of iteration (m+1) and using the 
constrained least square method, we search for the new solution k (m+1) that minimises the following function: 

|| x (m+1) Jfc (ffl+1) - b || 2 (8) 

submitted to the constraints of equation[6]and where b is defined by the relation b - W i, s - W sy „ + Xk (fl 'K 

The solution of the constrained least square method is obtained using procedures from the NAG library. 
3.3. The constraints: 

As the spectral resolution and the wavelength range are limited, the number of uncorrected stars is also limited. This leads to 
inherent incompleteness of the database which together with different photometric accuracy between galactic and stellar spectra 
can lead to non-physical solutions. So, one may have to constrain the solutions, in particular, in such a way that the Initial Mass 
function of stars (IMF) satisfies the known shape of this function derived from observations of resolved objects. 



Mass group 


Mass interval 


stars 


1 


17M G - 3OM 


O7-B0V 
M2Ia 


2 


3M Q - 17M Q 


B3-4V 

G0Iab,K4Iab,rG2Iab,rK0II,rK3Iab 


3 


1.6M -3M 


A1-3V 


4 


0.8M G - 1.6M Q 


F2V,F8-9V,G4V,rG0lV,rG5lV 
G0-4III,wG8III,G9III,K4III,M0.5III,M4III,M5III,rG9III,rK3III,rK3III(bis),rK5III 


5 


< 0.8M o 


G9-K0V,K5V,M2V,rK0V,rK3V,rMlV 



Table 1. The mass groups cutout of the stellar database in the case of the "Standard mode". 



We first define groups of main sequence (MS) stars following their range of lifetime. This translates into stellar mass 
intervals. In each of these groups, we include the evolved stars of the same initial mass. Tables and |2] show the cutout of the 
stellar database in "mass groups" for two different mass resolutions. The mass interval associated to each star of the base was 
determined using Schmidt-Kaler tables (1982). When necessary, we interpolated the available masses in the tables. Concerning 
the evolved stars, we used the evolutionary tracks of Padova and Geneva groups to associate an interval of initial masses for the 
stars. Each star thus represents an evolutionary stage in a given mass group. We aim in any case at not privileging one of the 
classical IMFs (as the functions of Salpeter 1955, Scalo 1986 and Kroupa et al. 1993). 

We will discuss here two different modes of constraints corresponding to the two different samplings in mass of the H-R 
diagram shown in tables^and|2] 

The "Standard mode": 



In this mode, we impose a general constraint on the number of born stars of the different mass groups: 
N t < Nj 



(9) 
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Mass group 


Mass interval 


stars 


1 


17M - 30M o 


O7-B0V 
M2Ia 


2 


3M - 17M Q 


B3-4V 

G0Iab,K4Iab,rG2Iab,rK0II,rK3Iab 


3 


1.6M Q -3M G 


A1-3V 


4 


1.4M Q - 1.6M 


F2V 

M4ni,M5in 


5 


\.\M B - 1.4M Q 


F8-9V,rG0IV,i'G5IV 
G0-4III,wG8III,G9III,K4III,M0.5III,rG9III,rK3III,rK3III(bis),rK5III 


6 


0.8M Q - 1.1M 


G4V 


7 


0.7M Q - 0.8M o 


G9-K0V,rK0V,rK3V 


8 


0.5M Q - 0.7M o 


K5V 


9 


< 0.5M o 


M2V,rMlV 



Table 2. The mass groups cutout of the stellar database in the case of the "Decreasing IMF mode". 



where Ni and Nj are respectively the numbers of born stars with masses located in the mass intervals ]Mi a ,Mib[ and ]Mj a ,Mjt,[ 
and where Mjt, < Mj a . Even though a galaxy is not described by a single burst of star formation, equation|9]is valid as well for all 
the Main-Sequence (MS) stars produced by a continuous star formation scenario or by multiple bursts of star formation. As stars 
of higher masses evolve faster than the ones of smaller masses, equation|5]is also valid for all stars (MS and evolved) belonging 
to the same mass groups. 

This translates into the following constraints set: 

N M si<N M sj (10) 

N Ti <N Tj (11) 

where Nms i an d Nt ; are respectively the number of MS stars and the total number of stars (MS and evolved stars) in a given 
mass interval. In tabled the mass intervals are chosen the narrowest possible such that eauationsllOlandfTTIare satisfied in the 
case of the classical IMFs. 

A consequence of introducing such constraints is that the resulting solutions will present higher synthetic distances than in 
the unconstrained solution, and therefore the fit will be less good. 

The "Decreasing IMF mode": 

The second mode of constraints called "the Decreasing IMF mode" states that all initial mass functions are decreasing 
functions. Consequently, for a burst of star formation, the total number of born stars Nn of a given mass group (Nn = Nms i + 
Nevoivedi) of a given group of stars corresponding to a mass interval ]M ia ,Mib[ with a mean mass Mi (Mi = MM+Mik) i s suc h that: 

AlogM,- * AV ; 
Alog Mj 

where Mj < M t . 

Alog Mj and Alog Mj are respectively the logarithmic lengths of each mass interval. The difference with respect to equation 
8 is that the number of stars in each group is weighted by the mass. 

As in the "Standard mode", equation^lis valid for the number of MS stars on the one hand and the total number (MS and 
evolved) stars on the other hand. Then, one can write: 

Alog Mi * Nms i 

N MSi < A1 ' (13) 
Alog Mj 

Alog Mj * N T i 

N Ti < * ' . ' (14) 
Alog Mj 



where the notations are the same as for the "Standard mode" 
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As equation dl 2l i depends on the mass intervals, the resulting mass cutout is more refined for the low mass stars (see tabled- 
Because of the better resolution of the mass cutout in this case, the number of constraints is increased compared to the "Standard 
mode". Consequently, the synthetic distance in this mode of constraints will be higher. 

The previous description of the adopted constraints shows that these take the form of large inequalities between the numbers 
of stars or equivalently, in the optical domain, between the stellar contributions to luminosity at the reference wavelength k Aa . 
The relation between the number of stars of class i and their contribution to luminosity is as follows: 

»< = : (15) 

where i /!(J g a i is the galactic (or the cluster) luminosity at the reference wavelength Aq and l Ao . is the luminosity of a star of class i 
at the same wavelength. 

As the constraints have the form of large inequalities (i.e. equalities are allowed) one may obtain solutions on the border of the 
domain of constraints (i.e. if the constraints are satisfied with equalities), the solution is probably over-constrained. In such a 
case, one has to conclude that the constraint set is not "appropriate" to the data. 



4. Results 

The effect of adding constraints to the inverse method has been tested on the globular cluster G170 located in M31 and three 
central regions of galaxies corresponding to three types of galaxy activity considered in Paper IV. 

As the aim of the present paper is to study the effects of constraining the solutions, we will only concentrate on a description 
and on the comparison of the different solutions. Therefore no complete physical discussion (especially in the case of the 
galaxies) will be made in the following sections since astrophysical issues about these objects have been plenty discussed in 
Paper IV (for a full discussion, we invite the reader to refer to the paper). 

In tables[3]to[6] we show the list of stars present in the database (column one), the unconstrained solution (only statisfying the 
positivity condition (|2j and equation Q) in column 2. In the two following columns, we show the mass intervals corresponding 
to the mass cutout for the "Decreasing IMF mode" and the constrained solutions of this mode. Columns 5 and 6 exhibit the 
mass intervals and the constrained solutions in the case of the "Standard mode". Finally, in some examples, we show particular 
solutions where we imposed a fixed value to one or more stellar contributions. 
The solutions list the stellar contributions, to luminosity at the reference wavelength Aq = 5450A. 



4. 1 . The stellar cluster G 1 70 

In order to validate our method, we compare the constrained and unconstrained solutions obtained for a single stellar population 
system, namely, the well studied globular cluster G170 located in the galaxy M31. The globular cluster spectrum is taken from 
Jablonkaetal(1992). 

The unconstrained solution (see table|3 shows that this cluster has a solar metallicity. The turnoff in G4V suggests an age of 
about 10 1() years and the contribution of dwarf stars to the luminosity of 57% shows that the luminosity in the optical is dominated 
by the MS stars. The metallicity is clearly solar (the contribution of stars of solar metallicity is about ~ 92%). The measured 
internal reddening E(B-V) is of ~ 0.05. This value of the reddening is determined as the correction to be applied to the observed 
spectrum to match the synthetic one. In this process the Galactic law is parametrized as in Howarth (1983). 

These results agree well with those found by Jablonka et al. (1992); as a matter of fact, the authors concluded, using mea- 
surements of equivalent widths of absorption lines, that this cluster has a solar metallicity as well as an age of 2 10 10 years. 

The "Decreasing IMF mode" presents a solution with acceptable synthetic distance in the sense that its value is at ~ 1 cr from 
the synthetic distance value of the unconstrained solution, but the contributions are slightly different from the latter: only ~ 38% 
of the optical luminosity is due to MS stars and ~ 20% of it is due to metallic stars. 

The constrained solution of the "Standard mode" is equal to the unconstrained one. This result shows that the latter is 
physically acceptable. 

The difference between the two solutions, "Decreasing IMF mode" and "Standard" or unconstrained mode is illustrated in 
Fig.l. 



4.2. The nucleus of the LINER NGC 4278 

As shown in table @] column 2, the unconstrained model suggests that the optical spectrum of this region is dominated by dwarf 
stars as well as a metal rich population. An internal reddening E(B-V) of ~ 0.02 is detected. 

The solution of the "Decreasing IMF mode" presents a much larger synthetic distance (with values exceeding the one of the 
unconstrained solution by more than 5cr). This shows that this mode is very constraining for the galactic region, a fact that is 
confirmed through the large discrepancy between the observed and the synthetic spectrum in Fig. [3 
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G170: Unconstrained solution 



o 



1.2 
1 

0.8 
0.6 
0.4 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



J i i i i I i i i_ 



5000 



1.2 
1 

0.8 
0.6 
0.4 



5000 



6000 7000 
X (A) 

G170: Decreasing IMF mode 



8000 



"| 1 1 1 1 1 1 1 r~ 




-| i i r 



J I L_ 



_l I I L_ 



LI I L_ 



6000 



7000 

x (A) 



8000 




J i i i i I i i i_ 



-| i i r 




J i i i_ 



Fig. 1. NGC G170: Synthetic and observed spectra respectively in dark and light lines for different solutions. 



The "Standard mode" provides a solution equivalent to the unconstrained one. It shows a small contribution of G4V stars 
badly determined suggesting an earlier location of the turnoff. This led us to impose a value to the contribution of this spectral 
class and to obtain an acceptable solution (shown in the last column of table 0} presenting an earlier turnoff and satisfying the 
constraints of this mode. 



4.3. The nucleus of the starburst galaxy NGC3310: 

According to the unconstrained solution, the MS stars and the metallic stars contribute, respectively, to ~ 94% and ~ 68% of the 
luminosity in the nucleus of NGC 3310 (table |5}, implying a stellar population dominated by MS and metallic stars. The best 
defined turnoff is situated in A1-3V but a turnoff in O7-B0V is possible as well since a small contribution of these stars is present 
but badly determined. The reddening is very high in agreement with the location of the turnoff that indicates an important event 
of star formation and consequently the presence of a big quantity of dust in the region. 

The solution of the "Decreasing IMF" mode presents a synthetic distance a little more than ~ lcr higher than the one of the 
unconstrained solution. It may therefore be acceptable (see also Fig. [3}- This solution distributes the non zero contributions to a 
larger number of dwarfs and confirm the high contribution to luminosity of dwarfs and metallic stars (resp. ~ 95% and ~ 49%) 
as well as the location of the turnoff situated in A 1-3 V (corresponding to an age of 200 million years for the last burst of star 
formation). In this solution, the constraint involving the O7-B0V stars is satisfied on the border of the domain of constraints (i.e. 
with equalities). This suggests that the "Decreasing IMF mode"is very restricting in this object; the absence of the hot 07-BOV 
stars might therefore not be real. This conclusion is supported by the presence of emission lines in the spectrum of the starburst 
galaxy which suggests an ongoing star formation occuring in the nucleus of this galaxy. 

The same scenario occurs in the "Standard mode" where the solution provides an acceptable synthetic distance. The best 
defined turnoff is situated in G5IV but small non zero contributions (not well defined) show that a turnoff at earlier types is 
possible as well. We show in the last column of table [5]the example of such a situation where we impose the contributions of 
B2-3V star in addition to the "Standard mode" constraints. In this solution, the synthetic distance is acceptable which leads us to 
the same conclusion as previously. 
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Unconstrained 


Mass interval 


Dec. IMF 


Mass interval 


Standard 


Star 


solution 


Dec. IMF mode 


mode 


Standard mode 


mode 


07-BOV 





17M - 3OM 





17M - 30M o 





B3-4V 





3M Q - 17M Q 


0.01 ±2 


3Af Q - 17M 





A1-3V 





1.6M Q - 3M Q 





1.6Af -3Af 





F2V 





0.8M G - 1.6M Q 





1.4M G - 1.6M Q 





F8-9V 





0.8M o - 1.6M Q 





l.lAfo- 1.4M e 





G4V 


22 ± 11 


0.8Af o - 1.6M 


7± 12 


O.8M - 1.1M Q 


22 ± 11 


G9-K0V 





< 0.8M o 





0.7Af o - 0.8Af o 





K5V 


27 ±7 


< 0.8M o 


11 ±5 


O.5M - 0.7M o 


27 ±7 


M2V 





< 0.8M G 


2± 1 


< 0.5M o 





rGOIV 





0.8M o - 1.6M Q 





l.lAfo - 1.4M Q 





rG5IV 





0.8M o - 1.6M Q 





l.lAfo - 1.4M Q 





rKOV 


8±6 


< 0.8M Q 


18 ±5.5 


0.7M G - 0.8M G 


8±6 


rK3V 





< 0.8M Q 





0.7Af o - 0.8Af o 





rMlV 





< 0.8M G 





< 0.5M o 





G0-4in 


31 ± 13 


0.8Af Q - 1.6M Q 


50 ±7 


1.1M Q - 1.4M Q 


31 ± 13 


wG8III 





0.8M o - 1.6M S 


0.01 ± 6 


1.1M - l-4Af Q 





G9III 





0.8M o - 1.6M Q 





l.lAfo - 1-4M 





K4III 





0.8Af o - 1.6M 





l.lAfo - l-4Af Q 





M0.5III 


11 ± 1 


0.8Af Q - 1.6Af s 


8± 1 


l.lAfo - 1.4Af Q 


11 ± 1 


M4III 





0.8Af Q - 1.6Af s 





1.4Af Q - 1.6Af Q 





M5III 


1 ±0.1 


0.8Af Q - 1.6Af Q 


1 ±0.1 


1.4Af Q - 1.6Af Q 


1 ±0.1 


rG9III 





0.8M o - 1.6M Q 





l.lAfo - l-4Af 





rK3III 





0.8M o - 1.6M Q 





l.lAfo - l-4Af 





rK3IIIbis 





0.8M G - 1.6M Q 





l.lAfo - 1.4M 





rK5III 





0.8Af o - 1.6M 


2±2 


l.lAfo - l-4Af 





GOIab 





3M Q - 17M Q 





3Af Q - 17Af G 





K4Iab 





3M - 17M Q 





3Af Q - 17Af G 





M2Ia 





17M Q - 3OM 


1 ± 1 


17M - 30Af o 





rG2Iab 





3M Q - 17M G 





3Af Q - 17Af Q 





rKOII 





3M Q - 17M G 





3Af - 17Af Q 





rK3Iab 





3M Q - 17Af Q 





3Af Q - 17Af Q 





D orx 


7.6 ±0.4 




8.0 ±0.35 




7.6 ±0.4 


E(B-V) 


0.05 




0.08 




0.05 



Table 3. Results of the spectral synthesis of the globular cluster G170. Each column displays the stellar contributions to lumi- 
nosity at = 5450A with their standard deviations for the various solutions (with different modes and types of constraints) as 
well as for the unconstrained solution. The over-abundant and under-abundant stars are respectively designated by an "r" and a 
"w" preceding the spectral type and luminosity class. D 2 is the synthetic distance (or the "mean" residual EW) with its standard 
deviation computed in the appendix and E(B-V) is the reddening that can be deduced as described in the text. 



4.4. The nucleus of the Seyfert2 galaxy NGC2110: 

The population in the unconstrained solution is dominated by dwarf stars and is moderately metallic (~ 68% of the optical 
luminosity is due to MS stars and ~ 52% is due to overabundant stars, see tableland Fig.|5|i. The best defined turnoff is situated 
in K0V but according to the previous discussion, an earlier turnoff is possible as well. 

The solution of the "Decreasing IMF" mode is acceptable since its synthetic distance is slightly higher than the one of the 
unconstrained solution and lies at less than lcr from this one. In this mode, the contribution of dwarf stars is enhanced (~ 78%) 
while overabundant stars contribute less to the visible luminosity (only ~ 38%). The best defined turnoff is situated in K3V but a 
turnoff in F2V is also possible, a fact that is confirmed in the particular solution where we imposed a contribution of ~ 3% to the 
star class F2V (see also Fig. [5}- 

In the "Standard mode", the solution presents a synthetic distance equal to the one of the unconstrained solution. The solution 
is very similar to the unconstrained one. 

5. Conclusion 

The ideal case for a spectral synthesis giving a synthetic distance equal to zero would be the case where the signal to noise 
ratio of the galactic and stellar spectra goes to infinity and where the stellar database is complete. In such a case, all stars with 
spectral types later than the spectral type at the turnoff position would have non zero contributions to luminosity. But in practice 



o 
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NGC 4278: Unconstrained solution 
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NGC 4278: Decreasing IMF mode 
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5000 



8000 
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Fig. 2. NGC 4278: Same as Fig.E 



all unconstrained solutions show many zero contributions; this is due to the finite signal to noise ratio of our spectra and to the 
limitation of the stellar database which itself is due to the finite spectral resolution. 

Therefore, constraining the stellar population would a priori reduce the number of zero contributions because of the additive 
information introduced in this process. But as can be seen in the previous results, no large improvement in eliminating the 
zero contributions has appeared. This is probably due to the not perfect adequation of the observational data. Actually as the 
constraints are expressed by large inequalities (i.e. equalities are allowed) optima are usually located on the border of the domain 
in which solutions are constrained. 



The stellar synthesis method with constraints presented in this paper has been applied to the 27 regions of galaxies studied 
in Paper IV. In general, all 27 regions present "Standard mode" solutions equal or very similar to the unconstrained solution. 
Moreover, all zero contributions in the unconstrained solutions remain null in the "Standard mode" or have small ill-defined 
values and all well- and ill- determined contributions remain respectively well- and ill- defined. This result shows that "Standard 
mode" solutions are generally included inside the error bars of the unconstrained solution and when they are not, their synthetic 
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Unconstrained 


Mass interval 


Dec. IMF 


Mass interval 


Standard 


Particular 


Star 


solution 


Dec. IMF mode 


mode 


Standard mode 


mode 


solution 


07-BOV 





17M Q - 30M o 





17M G - 30M G 








B3-4V 





3M e - 17M G 





3M G - 17M G 








A1-3V 





1.6M Q - 3M 





1.6M G - 3M G 








F2V 





0.8M G - 1.6M S 





1.4M G - 1.6M G 








F8-9V 





0.8M o - 1.6M S 





l.lAf - 1.4M G 








G4V 





0.8M o - 1.6M Q 


1 ±2 


0.8M o - 1.1M Q 


0.2 ± 1 


7 


G9-K0V 





< 0.8M G 





0.7M G - 0.8M G 








K5V 





< 0.8M o 


42 ±4 


0.5M o - 0.7M o 








M2V 





< 0.8M G 


4±2 


< 0.5M o 








rGOIV 





0.8M o - 1.6M S 





1.1M G - 1.4M G 








iG5IV 





0.8M o - 1.6M S 





l.lMo - 1.4M G 








rKOV 


38+ 14 


< 0.8M o 





0.7M G - 0.8M o 


37 ± 15 


26 ± 15 


rK3V 


38 ±8 


< 0.8M G 


35 ±3 


0.7M o - 0.8M o 


38 ±8 


40 ±5 


rMlV 


2±2 


< 0.8M o 


5±4 


< 0.5M o 


2.5 ±2 


3±2 


G0-4III 





0.8M G - 1.6M Q 





1.1M G - 1.4M G 








wG8III 





0.8M G - 1.6M Q 





UMo- 1.4M g 








G9III 





0.8M o - 1.6M S 





I.IMq - 1.4M G 








K4III 





0.8M G - 1.6M G 





1.1M Q - 1.4M G 








M0.5III 


5±2 


0.8M o - 1.6M Q 





l.lMo - 1.4M G 


5±2 


4± 1 


M4III 


3.5 ±0.1 


0.8M G - 1.6M 


3 ±0.2 


1.4M - 1.6M S 


3.5 ±0.1 


3 ±0.3 


M5III 





0.8M o - 1.6M 





1.4M G - 1.6M G 








1G9III 





0.8M o - 1.6M Q 





1.1M - 1. 4M B 








rK3HI 


3 ± 6 


0.8M o - 1.6M S 





1.1M - 1.4M G 


3±7 


4± 6 


rK3IIIbis 





0.8M G - 1.6M S 





1.1M G - 1.4M G 








rK5III 


10 ±4 


0.8M G - 1.6M 


10 ±2 


l.lMo - 1.4M G 


10 ±4 


11 ±4 


GOIab 





3M Q - 17M G 





3M G - 17M G 








K4Iab 





3M e - 17M G 





3M G - 17M G 








M2Ia 





17M Q - 30M o 





17M G - 30M G 








rG2Iab 





3M G - 17M 





3M G - 17M G 








rKOII 





3M Q - 17M Q 





3M G - 17M G 








rK3Iab 





3M e - 17M 





3M G - 17M G 








D- or x 1 


12.3 ±0.95 




17.1 ± 1.6 




12.3 ±0.9 


13.0 ± 1.2 


E(B-V) 


0.02 




0.00 R 




0.02 


0.00 



Table 4. The different synthesis solutions for the nucleus of NGC4278. Same notations as in table[3] The underlined contribution 
in the particular solution is imposed. 

In the reddening row, R means that the synthetic spectrum is redder than the observed spectrum. 



distances are at several cr from that of the unconstrained solution. 

In the "Decreasing IMF mode', the number of star classes contributing to the synthesis is often larger than that of the 
unconstrained solution and of the "Standard mode". This fact affects especially dwarf stars and is due to the sharper distribution 
of stars in the mass groups of the H-R diagram in this mode. For the same reason and because the number of constraints is larger, 
the synthetic distances here are in general larger than those of the previous mode and of the unconstrained solution. 

In both modes ("Standard" and "Decreasing IMF") some solutions satisfy their constraints on the border of the domain. This 
shows that in such cases constraints are somehow too strong and induce bias. However, these solutions provide some indications, 
thanks to the error bars, allowing one to find acceptable solutions that satisfy the desired conditions inside the domain of 
constraints (see previous examples). 

As a matter of fact, the solution of the least square problem is the one that minimizes the synthetic distance; this happens 
often on the border of the domain of constraints but the goal is not to find the optimal mathematical solution, rather a "realistic" 
or physical one next to the minimum. 

All previous results are very well confirmed in the case of the globular cluster G170. This is a very important point since 
this object is constituted of a single burst of star formation; consequently, any deviation of the behaviour of the resulting stellar 
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ngc 3310: Unconstrained solution 




NGC 3310: Decreasing IMF mode 



1.2 



~i 1 r 



t — i — i — i — i — i — r 



5000 



6000 7000 8000 

ngc 3310: Particular solution 



~i 1 r 





Fig.3. NGC 3310: Same as Fig.E 



population due to the inclusion of astrophysical constraints can clearly be detected in this object. 



This study has shown that the inverse method described in this paper and in Papers I, II and III is very stable against the 
inclusion of additional astrophysical constraints, and is, therefore, very reliable. 



However, constraining the solutions and using the information provided by the error analysis allows one to find similar 
solutions with younger bursts of star formation. Thus, it is crucial to perform tests such as in the previous section and to discuss 
the results, especially the different possible locations of the turnoffs, i.e. the age of the last burst of star formation. 
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Unconstrained 


Mass interval 


Dec. IMF 


Mass interval 


Standard 


Particular 


Etoile 


solution 


Dec. IMF mode 


mode 


Standard mode 


mode 


solution 


07-BOV 


4±5 


17M Q - 30M o 





17M G - 30M G 


4 ±33 





B3-4V 





3M e - 17M G 





3M G - 17M G 


0.1 ±44.5 


10 


A1-3V 


21 ±9 


1.6M -3M 


30 ± 12.5 


1.6M - 3M 


22 ±26 


15 ±9 


F2V 





0.8M G - 1.6M S 


1 ± 17 


1.4M G - 1.6M Q 








F8-9V 


6± 6 


0.8M o - 1.6M S 





l.lMa - 1.4M G 


6±9 


6.5 ±6 


G4V 





0.8M o - 1.6M Q 


6± 11 


0.8M G - 1.1M Q 








G9-K0V 





< 0.8M G 





0.7M G - 0.8M G 








K5V 





< 0.8M o 


11 ±4 


0.5M G - 0.7M o 








M2V 





< 0.8M G 


2± 1 


< 0.5M o 








rGOIV 





0.8M o - 1.6M S 





1.1M G - 1.4M G 








iG5IV 


30+ 10.5 


0.8M o - 1.6M S 


27 ±9 


l.lMo - 1.4M G 


30 ± 10 


29 ± 10 


rKOV 


33 ±5.5 


< 0.8M o 


18 ± 7 


0.7M o - 0.8M o 


33 ±4.5 


34 ±5 


rK3V 





< 0.8M G 





0.7M G - 0.8M G 








rMlV 





< 0.8M o 





< 0.5M G 








G0-4III 





O.8M - 1.6M Q 





l.lMo - 1.4M G 








wG8III 





0.8M G - 1.6M Q 





l.lMo - lAM e 








G9III 





0.8M o - 1.6M S 





UMq - 1.4M G 








K4III 





0.8M G - 1.6M G 





1.1M Q - 1.4M G 








M0.5III 





0.8M G - 1.6M G 





l.lMo - 1.4M G 








M4III 


1 ±0.2 


0.8M o - 1.6M 


1 ±0.2 


1.4M G - 1.6M G 


1 ±0.2 


1 ±0.2 


M5III 





0.8M o - 1.6M 





1.4M G - 1.6M G 








1G9III 





0.8M o - 1.6M Q 





1.1M - 1. 4M B 








rK3HI 





0.8M o - 1.6M S 





1.1M - 1.4M G 








rK3IIIbis 





0.8M o - 1.6M S 





1.1M G - 1.4M G 








rK5III 


5±2 


0.8M o - 1.6M 


4±2 


l.lMo - 1.4M 


5±2 


5±2 


GOIab 





3M Q - 17M G 





3M - 17M G 








K4Iab 





3M e - 17M G 





3M G - 17M G 








M2Ia 





17M Q - 30M o 





17M G - 30M G 








rG2Iab 





3M G - 17M 





3M G - 17M G 








rKOII 





3M Q - 17M 





3M G - 17M G 








rK3Iab 





3M e - 17M 





3M G - 17M G 








D- or x 1 


11.4 ± 1.5 




13.2 ± 1.5 




11.4 ± 1.2 


11.6 ±1.5 


E(B-V) 


0.23 




0.2 




0.23 


0.23 



Table 5. Results of the spectral synthesis of NGC 3310 nucleus. Same notations as in tables[3]and|3] 



Appendix A: Calculus of the standard deviation on the synthetic distance: 

In this appendix, we compute the standard deviations on the synthetic distance due to observational errors around the studied 
object. This calculation is complementary to the error analysis made in Paper III where only the standard deviations around the 
stellar contributions and the variance-covariance matrices were computed. As the synthetic distance is a scalar, its variance- 
covariance matrix is reduced to its variance. 

Thus, here we search for the deviations d(D 2 ) around the square of the obtained synthetic distance due to deviations 
dW bs around the observation W o bs0- We recall that this computation is only valuable in the overdetermined case (in the 
underdetermined case, this distance is equal to zero). 

If we make a change of variables on the equivalent widths as W'j = P l / 2 Wj, the square of the synthetic distance will be written 
as follows: 

D 2 = (W' m - W' obs f(W' m - W' obs ) (A. 1) 

Then a differenciation around D 2 Q gives: 

dD 2 = 2(dW' syn - dW' oh f(W' syn0 - W'obso) • (A.2) 
Now, if we replace dW' syn by HdW' gbs , where H is the orthogonal projector on the synthetic surface (see Paper III), then we get: 

dD 2 = 2((H - ld)dW' ob f(W' syn0 - W' obs0 ) fA ~ 

= 2dW' T ohs (H T -\d)(W' syn0 -W' obs0 ). ( ^' 
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J I I L 
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NGC 2110: Decreasing IMF mode 
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5000 



8000 



Fig. 4. NGC 2110: Same as Fig.E 



where Id is the identity matrix. 

Let us write, on the one hand, the definition of the variance of the quantity dD 2 : 



varidD 2 ) =< dD 2 dD 2T > - < dD 2 > 2 



(A.4) 



on the other hand, we have < dD 2 >= 2Do < dD >. Then if we translate the origin of the vector space in W svn o, we can 
consider the subspace of dimension 1 of which the generator vector is W' o bs0 ~ W' svn Q. In this subspace the synthetic distance is 
described by the same vector Do = W' b s o - W synO and the deviation to this distance is vector dD — D - Dq. Thus, if we call u 

W W 

'* oiso rr „„o \ „,„ +, + „,.fU„„„„„i „„ n _ ,.,,T . 



the unit vector of this subspace (u - 



\\W' ab!0 -W' { 



-), we can construct the orthogonal projector over it as P = uu ; then we get 



dD = PdW obs and < dD >=< dD >= P < dW' sy „ >= 0. This implies that < dD 1 >= and var(dD 2 ) = var(D 2 ). 
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NGC 2110: Particular solution 



1.2 



1 

1 - 





o 

1/3 



0.6 - 



0.4 - 



5000 



6000 



7000 



8000 



X (A) 



Fig. 5. NGC 2110: Particular solution. Same as Fig. [I] 

Back to equation JA.4> . we have: 

var(dD 2 ) = var(D 2 ) 

=< dD 2 dD 2 > 

= 4AW' T (H - Id) < dW' obs dW' T obs > (H T - \d)AW 
where AW = W' obs0 - W' m o- Then 



In addition, as < dD >=< dD 2 >- 0, we set the simple equation 
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Unconstrained 


Mass interval 


Dec. IMF 


Mass interval 


Standard 


Particular 


Etoile 


solution 


Dec. IMF mode 


mode 


Standard mode 


mode 


solution 


07-BOV 


3±7 


17M Q - 30M o 





17M G - 30M G 


2 ±67 





B3-4V 





3M e - 17M G 





3M G - 17M G 


0.1 ±71 





A1-3V 





1.6M Q - 3M 





1.6M G - 3M G 








F2V 





0.8M G - 1.6M S 


0.006 ± 7 


1.4M G - 1.6M Q 





3 


F8-9V 





0.8M o - 1.6M S 





1.1 M Q - 1.4M G 








G4V 





0.8M o - 1.6M Q 


3± 18 


0.8M o - 1.1M Q 





2 ±14 


G9-K0V 





< 0.8M G 





0.7M G - 0.8M G 








K5V 


9±8 


< 0.8M o 


32 ± 6.5 


0.5M G - 0.7M G 


9 ±10 


32 ±7 


M2V 


6±2 


< 0.8M G 


6± 1 


< 0.5M o 


6± 1 


6± 1 


rGOIV 





0.8M o - 1.6M S 





1.1M G - 1.4M G 








iG5IV 


9± 12 


0.8M o - 1.6M S 


10 ± 16 


UMq - 1.4M G 


10 ± 12 


9+15 


rKOV 


19 ± 6 


< 0.8M o 





0.7M G - 0.8M G 


19 ±8 





rK3V 


22 ±5 


< 0.8M G 


27 ±3 


0.7M o - 0.8M o 


22 ±7 


26 ±2 


rMlV 





< 0.8M o 





< 0.5M G 








G0-4III 





O.8M - 1.6M Q 





l.lMo- 1.4M G 








wG8III 





0.8M o - 1.6M Q 





UMo- 1.4M g 








G9III 





0.8M o - 1.6M S 





UMq - 1.4M G 








K4III 


15 ±5 


0.8M o - 1.6M S 


6±4 


1.1M - 1.4M Q 


15 ± 7 


6±4 


M0.5III 


12 ± 3 


0.8M o - 1.6M 


12 ±3 


l.lMo - 1.4M G 


12 ± 1.5 


11.5 ±3 


M4III 


2 ±0.3 


0.8M o - 1.6M 


2 ±0.2 


1.4M - 1.6M G 


2 ±0.1 


2 ±0.2 


M5III 





0.8M o - 1.6M 





1.4M G - 1.6M G 








1G9III 





0.8M o - 1.6M Q 





1.1M - 1. 4M B 








rK3HI 





0.8M o - 1.6M S 





1.1M - 1.4M G 








rK3IIIbis 





0.8M o - 1.6M S 





1.1M G - 1.4M G 








rK5III 


2.5 ±7 


0.8M o - 1.6M 


1±8 


l.lMo - 1.4M 


2.5 ±4 


2± 8 


GOIab 





3M Q - 17M G 





3M - 17M G 








K4Iab 





3M e - 17M G 





3M G - 17M G 








M2Ia 





17M Q - 30M o 





17M G - 30M G 








rG2Iab 





3M G - 17M 





3M G - 17M G 








rKOII 





3M G - 17M 





3M G - 17M G 








rK3Iab 





3M e - 17M 





3M G - 17M G 








D- or x 1 


13.1 ±2.0 




13.5 ± 1.9 




13.1 ± 1.2 


13.5 ± 1.9 


E(B-V) 


0.05 




0.0 




0.05 


0.05 



Table 6. Results of the spectral synthesis of NGC 2110 nucleus. Same notations as in tables[3]and|3] 
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